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ABSTRACT 



This study presents a simplified approach to studying the dynamics of global 
texture collapse. We derive equations of motion for a spherically symmetric field 
configuration using a two parameter ansatz. Then we analyse the effective potential 
for the resulting theory to understand possible trajectories of the field configuration 
in the parameter space of the ansatz. Numerical results are given for critical wind- 
ing and collapse time in spatially flat non-expanding, and flat expanding universes. 
In addition, the open non-expanding and open-expanding cases are studied. 



1. Introduction 



Global texture theory is a cosmological model of large-scale structure forma- 
tion. Texture is a semi-topological defect in a theory where a global symmetry 
group G is broken to a group H such that ^3(7^) 7^ 1- We can take as a toy model 
a theory with a complex doublet of scalar fields and a Mexican hat potential gen- 
eralized to four dimensions, in which case the vacuum manifold is S s . 

Texture is formed when the field wraps around enough of S 3 to cause unwind- 
ing. Energy lumping of the field configuration causes accretion of matter, thus 
forming large-scale structure 111 121 . 

To understand cosmological structure formation, we need to have information 
about the texture distribution. The texture distribution depends crucially on the 
critical winding 131 (i.e. the winding above which collapse occurs). The formation 
probability for texture with a given winding distribution depends on the critical 
winding. More or fewer textures are formed if the critical winding is low or high, 
respectively. 

Critical winding is determined by the dynamics of the field configuration as 
uncorrelated regions of the field come into the horizon. Thus, we need to have 
a good understanding of field dynamics in the texture model to understand the 
texture distribution. 

Dynamics of texture field configurations has been analyzed in a variety of 
studies. A self-similar a-model solution has been found for the flat, non-expanding 
universe with winding w — l m . For this solution collapse proceeds at the speed of 
light. This study is limited by the fact that dynamics at the point of unwinding 
are not obtainable with the approximations used. Also, solutions for non-integer 
winding are not known. 

Numerical studies have been able to get past the above limitations 141 [5) 161 . The 
full equations have been integrated for spherically symmetric and ellipsoidal config- 
urations in flat, expanding universes; non-integer winding and the dynamics at the 
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point of collapse were also included in these studies. Furthermore, full cosmological 
simulations have also been done, including baryonic matter effects . 

Some analytical work has been done in the adiabatic approximation using the 
cr-model approach . These studies have approximated the Hamiltonian for large 
R, where R is a physical cutoff understood to be half the intertexture separation, 
and used a two parameter ansatz for a spherically symmetric field configuration. 

This study extends the above analytical work to include kinetic energy. Here 
we investigate a constrained Hamiltonian, where the field configuration is (as in 
the above study) approximated by a two parameter ansatz. 

First, we investigate the shape of the effective potential due to gradient energy 
to understand what sort of trajectories can be expected. Then, we integrate the 
Hamiltonian equations of motion to find the trajectory followed by the field ansatz 
in phase space. 

Analysis of the effective potential for this theory results in a clear understanding 
of why trajectories in phase space evolve as they do. A saddlepoint is found lying on 
the barrier between expanding and collapsing trajectories explaining the existence 
of trajectories which appear initially to be headed for collapse, but then expand, 
and trajectories which appear to be headed toward expansion, yet finally collapse 
(see section 3). 

Critical windings were found for fields uncorrelated on scales larger than the 
horizon for intertexture separation 2R, where R is in units of the horizon size, to 
be: 

For R = 2.0 (For it! = f .5) (errors are plus or minus the difference to the next 
closest calculated trajectory.) 

0.6562 ±0.0001 (0.6967 ±0.0001) (flat case), 
0.6693 ± 0.0001 (0.7471 ± 0.0001) (radiation era), 
0.667 ± 0.007 (0.7543 ± 0.0001) (matter era). 

Collapse time was found to increase at critical winding when approached from 
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high winding. This phenomenon is also explained by the saddlepoint found in the 
analysis of the effective potential. In fact, critical winding should, according to 
arguments based on the shape of the effective potential, go to infinity at critical 
winding. 



2. Derivation of Equations of Motion 

As was mentioned in section 1, global texture can be described by a theory 
with a complex doublet of scalar fields with a Mexican hat potential. The action 
for this theory is 



S = J (fx[dv§ a d^ a -\(\§\ 2 -rf)}^ 

where a = 1, 4 and rj is the scale of symmetry breaking. 

For low temperatures relative to the phase transition and times before un- 
winding we can fix the scalar field at the minimum of the potential and treat the 
potential as a constraint. This gives us a a-model action 

S = J d 4 x[d^^ a d^ a }^ 

with the constraint 

l$| 2 = ri 2 - 



We can use the spherically symmetric ansatz 



$ a = i](smx sin 6* sin 0, sin x sin 9 cos0, sin x cos 9, cosx), 



where x — xij'i t) is a radial field and the background Friedmann- Robertson- 
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Walker metric 



a 



2 



9flu = diag(l, 



1 — kr 2 



2 2 2 2 • 2 n\ 
, —a r , —a r sin 



to obtain the action 



S 




, 4 r , 2 x /2 (l-^ 2 ) 2sin 2 x 1 a 3 r 2 



a 2 r 2 J v /T3^ 



In the flat (A; = 0) non-expanding case (a(t) = 1), the solution to the equations 
of motion is 



and vary £ to obtain non-integer winding. In this ansatz a = a(t) and £ = £(£), 
£ is the point of intersection of the two segments of the ansatz and R is a cutoff, x 
is initially taken to be the horizon size (the scale on which the field is correlated) 
and R is taken to be 1/2 the intertexture separation, a fixed multiple of if -1 . The 
initial horizon H(to) is set to be equal to 1. 

This ansatz is expected to be good for understanding precollapse dynamics, 
especially around the time when the field configuration enters the horizon. 

With this ansatz in the above action, we derive the Hamiltonian 
H = p^aV/i)- 1 +p|(4a 3 7r 2 / 4 )" 1 

+aa 2 ir 2 I 2 + 2a/ 3 + 2a sin 2 (£7r)/ 5 

where p a and are momenta canonical to a and £ respectively, and 1% are integrals 
over the radial variable r. 



X = 2arctan(— |) for t < 0. 
We can mimic the exact solution with the simple ansatz (Figure 1) 



X = caxr for < r < x 
X = for x < r < R 
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From this Hamiltonian we find the equations of motion in phase-space (p a , p^, a, £)| 
to be 



d = p a (2a 3 n 2 I 1 )- 1 
i = p^2a\ 2 h)- 1 



-2ao;7r i2 — 7r — 2a— — 

oa oa 



-2asin 2 (^)|| 



p, = pi (4aV)- 1 / r 2 ^ + p| (4aV)- 1 / r 2 ^ 



2 2 <9/ 2 9/ 3 



-47ra sin(^7r) cos(^7r)/5 — 2a sin (^ 7r )"^' 



and the integrals /j are 

/i = / dr 



r 4 







\f\ — kr 2 



X 

h = J drr 2 y/ 1 — kr 2 



X 

f , sin 2 (a7rr) 
73 = / dr ; 





^,2 

I4 — I dr 



\/l — kr 2 



R 

dr 



yl — kr 2 



X 



The integrals are all trivial except ^3 in the open case (k ^ 0). For this case, 
we expand sin 2 r to sufficient order and calculate the resulting integrals. To solve 
the equations of motion in the open case, we need expressions for a(t), the scale 
factor of the universe. We use the fact that the universe is close to flat up to its 
current point of evolution and use the flat space radiation- and matter-dominated 
expanding universe scale factors as approximations to the actual scale factor. 

Finally, the expression relating £ to winding w is 

sin(27rO 



w(0 = £ 



2tt 



3. Results 

The terms in H independent of momenta p a and form the effective potential 
due to gradient energy in the field. Here (Figure 2) we plot a portion of it for the 
flat non-expanding case with R — 1.5. 

There are two attractors in the potential, one in the upper righthand corner 
(large a, large £) to which collapsing configurations approach; and one in the lower 
lefthand corner (small a, small £) to which expanding configurations approach. 

Notice the saddlepoint. To the lower left (small a, small £) and upper right 
(large a, large £) are the attractors for collapsing and expanding configurations. 
To the upper left (small a, large £) and lower right (large a, small £) we have 
regions of high potential. 
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If R is increased, the saddlepoint moves leftward to small values of a with £ 
remaining approximately the same, and its size shrinks (that is the gradient gets 
steeper on all sides). If R is decreased, the saddlepoint moves to the right to 
larger values of a with £ remaining approximately the same, and its size increases 
(gradients get less steep on all sides). The gradients around the saddlepoint get 
large for small R because the total energy of the configuration increases as R gets 
large and thus the height of the barrier increases. 

We take initial conditions for £ and a based on the physical constraint x^niai = 
1. This means that the initial £ and a obey the relation £ = a. For intertexture 
separation 2R = 4, this implies that for winding near critical the starting point of 
the field evolution trajectory is on the upper left region (hill) of the saddlepoint. 

For our initial conditions, x initial = 1 and thus a ~ p a and £ ~ p^, so near 
the saddlepoint, the field momenta can be thought of as 'rolling' on the effective 
potential (farther away momentum behavior becomes more complex). 

Therefore the field will roll down the hill and up the barrier on the opposite 
side of the saddle. Depending on whether the initial configuration is above or 
below critical winding the trajectory will scatter off the upper or lower side of the 
righthand barrier. If the trajectory moves to the basin in the upper right of the 
potential, the configuration collapses and unwinds. But if the trajectory scatters off 
the lower part of the righthand barrier, the configuration collapses for some time, 
but then, as the configuration falls to the lower left attractor, the configuration 
expands. 

For large R, the opposite effect occurs. Initially trajectories start on the lower 
right hill, then they roll across toward the opposite barrier and scatter to one side 
or the other depending on initial winding. These trajectories were discovered by 

rgi 

Perivolaropoulos in the large R limit. 

Trajectories that don't start on the saddlepoint fall more or less directly into 
the attractors on their respective sides of the saddlepoint. 
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For large R, there are fewer trajectories that start on or near the saddlepoint 
because its size shrinks. 

[3] 

Monte Carlo simulations suggest that , if we assume random correlations on 
length scales larger than the horizon, intertexture separation should be around two 
times the horizon size. This indicates that, usually, few textures will expand then 
collapse, but many will first tend towards collapse then expand, due to the location 
of the saddlepoint with respect to typical initial conditions. 

The following figures show integrated trajectories with Xi n m a i = 1 for flat 
non-expanding, flat expanding and open expanding universes. The trajectories are 
plotted in (£, a)-space. 

Critical winding can be found either by distinguishing which trajectories finally 
expand and which finally collapse, or by calculating the intersection of the line of 
initial conditions and the line of intersection of the tangent surface to the top of 
the ridge in which the saddlepoint lies with the ridge. We found the former method 
more convenient. 

Note that critical winding increases from the non-expanding (Figure 3) to the 
radiation-dominated (Figure 4) to the matter-dominated cases (Figure 5). This 
can be understood physically. The expanding background introduces an extra 
pull on the texture configuration, thus causing more configurations to expand, 
and thus pushing critical winding higher than in the non-expanding universe. The 
matter-dominated background expands more quickly than the radiation-dominated 
background, thus pushing critical winding even higher in the matter-dominated 
case. 

The effect of extra pull due to the background manifests itself in the effective 
potential such that the barrier region, including the saddlepoint moves to a higher 
value of £. 

In an open background, the area of concentric spheres increases more quickly 
as we leave the origin than in a flat background. This effect also tends to add a 
pull to the configuration, and push critical winding upwards. 
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As R increases (Figure 6), the saddlepoint in the effective potential moves to 
the left in the figures. This causes the line of initial conditions to move across 
the saddlepoint and critical winding (the point of intersection of the line of initial 
conditions and the top of the barrier) goes down, since the saddlepoint is in a 
skewed orientation with the lefthand region of high potential at larger £ than the 
righthand region of high potential. 

Critical windings were found for fields uncorrelated on scales larger than the 
horizon (x initial = 1) for intertexture separation 2R, where R is in units of horizon 
size, to be: 

For R = 2.0 (for R = 1.5) (errors are plus or minus the difference to the next 
closest calculated trajectory.) 

0.6562 ±0.0001 (0.6967 ±0.0001) (flat case), 
0.6693 ± 0.0001 (0.7471 ± 0.0001) (radiation era), 
0.667 ± 0.007 (0.7543 ± 0.0001) (matter era). 

Collapse time is affected by the length of time the configuration takes to roll 
off the saddle. Configurations very near to the saddlepoint take longer to roll off, 
since the saddlepoint is a point of unstable equilibrium. For configurations starting 
at the unstable equilibrium point, collapse time should be infinite. 

Also, there are trajectories below critical winding that evolve for a substantial 
period of time as if they would collapse, but eventually expand. Although there 
is no definable collapse time for these configurations, there might be appreciable 
matter accretion possible. 

Collapse time, defined as the time it takes for x, the joining point of the two 
segments in the ansatz, to reach half its initial value, for the flat non-expanding 
universe is plotted in Figure 7. 
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4. Discussion and Concluding Remarks 



We have used a two parameter spherically symmetric ansatz in the low temper- 
ature a-model approximation to the texture action to obtain a Hamiltonian and 
equations of motion for texture dynamics. 

From the Hamiltonian we have isolated the effective potential and found a sad- 
dlepoint that explained the behavior of the field configuration. The saddlepoint 
explains configurations which expand then collapse, and configurations which col- 
lapse then expand, as well as simpler trajectories in which the field configuration 
simply collapses or simply expands. 

From the Hamiltonian we derived equations of motion which were ordinary 
differential equations, in contrast to partial differential equations. The ODEs were 
easy to integrate on the computer, and results were obtained ~ 1000 times faster 
than integrating the PDEs. 

Critical windings were found for fields uncorrelated on scales larger than the 
horizon {xi n m a \ = 1) for intertexture separation 2R, where R is in units of horizon 
size, to be: 

For R = 2.0 (For R — 1.5) (errors are plus or minus the difference to the next 
closest calculated trajectory.) 

0.6562 ±0.0001 (0.6967 ±0.0001) (flat case), 
0.6693 ± 0.0001 (0.7471 ± 0.0001) (radiation era), 
0.667 ± 0.007 (0.7543 ± 0.0001) (matter era). 

Collapse time was found to increase as winding approached critical from above 
(for collapsing configurations). 

This study is also interesting, in general, as an approach to studying the a- 
model for non-integer windings where there is a cutoff scale R. 
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FIGURE CAPTIONS 



1) The ansatz for x as a function of r 

2) A contour plot of the effective potential. 

3) Trajectories are projected on (a, £)-space. In this figure trajectories are plot- 
ted for the spatially flat non-expanding universe with R = 1.5, x = 1.0. All 
trajectories start on the line £ = a 

4) Trajectories are plotted for the radiation-dominated universe with R — 1.5 
and x = 1.0. 

5) Trajectories are plotted for the matter-dominated universe with R— 1.5 and 
x = 1.0. 

6) Trajectories are plotted for the flat case with R = 2.0 and x = 1.0. 

7) Collapse time, defined as the time it takes for x to reach half its initial value, 
is plotted versus winding, for a non-expanding universe. 
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